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Abstract. Type la supernovae (SNe la) are a prime tool in observational 
cosmology A relation between their peak luminosities and the shapes of their 
light curves allows to infer their intrinsic luminosities and to use them as distance 
indicators. This relation has been established empirically However, a theoretical 
understanding is necessary in order to get a handle on the systematics in SN la 
cosmology Here, a model reproducing the observed diversity of normal SNe la 
is presented. The challenge in the numerical implementation arises from the vast 
range of scales involved in the physical mechanism. Simulating the supernova on 
scales of the exploding white dwarf requires specific models of the microphysics 
involved in the thermonuclear combustion process. Such techniques are discussed 
and results of simulations are presented. 



1. Introduction 

SNe la are extremely bright cosmic explosions with properties that are more 
homogeneous than those of other astronomical transients. Furthermore, a cor- 
relation between the width and the decline rate of the -B-ba nd light curve 
( "width-luminosity relation", WLR) allows to calibrate them ( Phillips! 1 19931 ; 



Phillips et al. 19991 ) and makes them the best distance indicators out to red 



shifts of about one. This correlation, however, is established only empirically 
on a set of nearby SNe la. Consequently, a quantification of systematic errors 
resulting from the calibration procedure is difficult to achieve. Moreover, aiming 
at distance determinations to far events, evolutionary effects could potentially 
obscure the measurements. It is clear that a sound understanding of the physics 
of SNe la is desirable in order to improve their precision as distance indicators 
in observational cosmology. 

Besides adding to the ever growing cosmological SN la databases, obser- 
vations in the past decade allowed to take a close look at a number of nearby 
events. It turned out that, although they form a remarkably homogeneous class , 



Mazzali et al. 


2007 


Branch et al. 


1993) 



. Apart from variations within the "normal" (as defined by 
SNe la, there are distinct sub-classes which differ signifi- 
cantly from the bulk of events. A first step towards a physical understanding of 
the class of SNe la and to improve their quality as distance indicators is therefore 
to identify the origin of this diversity. 
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A model that reproduces a large range of observational characteristics will 
be discussed in the following. Recent numerical simulations suggest that it 
potentially can account for the "normal" SNe la. This, however, also implies 
that the distinct sub-classes have to be explained in different physical scenarios. 



2. Astrophysical Model 



SNe la are attributed to thermonuclear explosions of carbon-oxygen white dwarf 
(WD) stars. In order to evolve a WD to a state where a thermonuclear explosion 
can trigger, the supernova progenitor has to be a binary syst em. Different scenar- 



ios of the progenitor evolution have been suggested (see, e.g.. lHillebrandt Sz Niemever 
2000l ) and may account for different SN la subclasses. Here, we focus on the 
so-called Chandrasekhar-mass explosion scenario, where the WD has accreted 
matter from its companion so that its mass approaches the stability limit - 
the Chandrasekhar mass of ~1.4M Q . Close to this mass limit, the density 
in the core of the WD increases dramatically so that carbon burning ignites. 
This, however, does not yet trigger the actual explosion process because convec- 
tive cooling still moderates the b urning. The resulting pre-explosion simmering 
phase lasts for about a century ( Wooslev et al.l [20041 ) and is characterized by 
highly turbulent convective motions. Gradually, the background temperature 
increases so that finally one or many hotspots in the turbulent flow near the 
WD's center undergo a thermonuclear runaway. Out of these "ignition sparks", 
a flame starts to propagate. It incinerates the WD material and leads to an 
explosion of the star on time scales of 1-2 s. As discussed in Sect. 23 the flame 
propagation is highly sensitive to the initial conditions and the ignition geome- 
try strongly affe cts the strength of the overall explosion process. Despite some 
recent progress (Hoflich fc Stein 20021; Wooslev et al.ll2004l ; IWunsch fc Wooslev 
20041 ; iKuhlen et al.ll2006l : IZing ale et al.ll2QQ9l i7 the physics of the ignition process 
and the turbulent simmering phase is still extremely challenging and realistic 
conditions are out of reach for numerical implementations. Therefore, the num- 
ber and spatiotemporal distribution of ignition sparks remains uncertain. 

After ignition near the center of the WD, the thermonuclear burning front 
propagates outwards. Two distinct mo des of propagation , a su bsonic deflagra- 
tion and a supersonic detonation (see Landau fc Lifshit3 19591 ) are consistent 
with the conservatio n laws of hydrodynamics. The first attempts to model- 
ing SN la explosions (Arnett 19691 ; lArnett et ailll97ll ). however, showed that a 
prompt detonation converts the WD almost entirely into iron-group elements. 
This is in contradiction with the observed spectra of these events which show 
strong intermediate-mass element (IME) features. Such elements (like Si, Ca, 
and S) are synthesized in thermonuclear explosive burning at low densities and 
therefore the WD material must expand prior to incineration. This can only 
be achieved if the flame propagates subsonically, i.e. as a deflagration. Before 
reaching low densities, however, the flame burns the high-density core material 
to iron-group elements, predominantly 56 Ni. In its radioactive decay, this iso- 
tope releases Gamma-rays which are scattered down to optical wavelengths in 
the ejecta and make the supernova bright. 

A laminar deflagration is a very slow process (see Timmes <fe Wooslev 19921 . 
for laminar flame speeds). In order to burn sufficient amounts of material to 
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explode the WD, the flame propagation must accelerate significantly. And in- 
deed, such an acceleration is to be expected as the flame propagation from 
the WD's center outwards produces a buoyancy-unstable stratification of light 
and hot ashes under dense fuel in the gravitational field of the star. The en- 
suing Rayleigh- Taylor instability leads to shear flows between rising plumes of 
burning material and fuel downdrafts with Reynolds numbers of ~10 14 . As a 
consequence of shear instabilities, a turbulent eddy cascade establishes. Down 
to the Gibson scale at which the laminar flame speed equals the eddy velocity, 
the flame interacts with turbulent eddies of this cascade. For most parts of the 
explosion, the Gibson scale lies orders of magnitude above the flame thickness 
which is only of the order of millimeters to centimeters. Therefore, turbulent 
eddies wrinkle and corrugate the flame on large scales without affecting the mi- 
crophysics of the burning. In this flamelet regime of turbulent combustion, a 
sufficient flame acceleration is achieved to explode the WD (e.g. 



Reinecke et al.l 



l2002l : iGamezo et a"Dl2003l : IRopke k Hillebrandtll2005l : IRopke et alJl2007ah . Only 
when the fuel density falls below ~10 7 gem -2 (due to the WD expansion and in 
the outer layers of the star), the Gibson scale becomes smaller than the broad- 
ening flame structure. In this late stage of the explosion, turbulent edd i es sti r 
the internal flame structure mixing fuel and ash. Niemever Sz Hillebrandt ( 19971 ) 
suggested that this may cause the transition from an initial subsonic d eflagration 
to a supersonic detonation in t he delayed detonatio n SN la scenario (iKhokhlovi 
Il99ll ). Microphysical studies (jLisewski etail |2000| ; fWooslev et ail 120091 ) iiidT 
cate that entering the distributed burning regime alone may not be sufficient 
for causing a deflagration-to-detonation transition (DDT), but in addition tur- 
bulent velocities of ~1000kms~ 1 at scales of about 10 km are required in this 
late bur ning regime. Analyzing three-dimensional simulations of deflagrations 
in WDs, Ropke ( 20071 ) found that such high turbulent fluctuations do occur in 
late stages of the explosion. 



3. Numerical Techniques 

The challenge in the numerical implementation of the above astrophysical sce- 
nario of thermonucl ear supernova explosions ar i ses from the multi-scale character 
of the problem (e.g. lRopke Bruckschenir2008l ; lRopkell2008l ). The physically rel- 



evant range in spatial scales covers about 11 orders of magnitude. In particular, 
the representation of a thin thermonuclear flame on a grid comprising the entire 
WD and the flame / turbulence interaction r equire special numerical approaches 
(for a summary see IRopke fc Schmidtll2009l ). 



In the implementa tion used here, the fla me propagation is modeled based on 
the level-set technique ( Reinecke et al. 19991 ). It is associated with the zero-level 



set of a signed distance function G which is defined positive in the ashes and neg 
ative in the fuel. The advancement of the flame is described by a partial differ- 
entia l equation modifying this G- field in an appropriate way (see iReinecke et al 



19991 ). As an input, this requires the flame propagation speed, which, in the 



flamelet regime, is set by the turbulent velocity fluctuations on the scale of 
computational grid. These are derived from a sub grid-scale turbulence model 
(|Niemever & Hillebrandtlll995l : ISchmidt et alll2006l ^ 
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Following lGolombek &: Niemeverl ( 20051 ) and lRopke k, Niemeverl ( 20071 ). the 
detonation wave is propagated by a separate level-set function. This allows to 
prevent the detonation fro m unphysically crossing ashes left behind from the 
previous deflagration stage ( Maier fc Niemever 20061 ). 

Another challenge arises from the expansion of the WD and the SN la 
ejecta. In order to reliably compute synthetic observables from the results of 
explosion simulations, these have to be followed to a state of hydrodynamically 
relaxed homologous expansion. This requires to simulate the explosion process 
for at least 10 s. An comp utational gri d co-expanding with the ejecta allows to 
keep then on the domain ( Ropkel [20051 ) and with moving nested gr ids the flame 



front can be optimally resolved for a fixed number of grid cells (jRopke et al 
1200631 ). 

Another scale problem arises from the discrepancy between hydrodynamic 
time scales and the time scales of some of the involved nuclear reactions. Cur- 
rently, this problem is bypassed in our implementation by only following a very 
limited set of nuclear species in the explosion simulation a nd reconstructing the 
details of the nucleo synthesis in a post-processing step (e.g. lTravag lio et alJl200l 
Ropke et al.ll2Q06al ) . This is achieved on the basis of passive tracer particles ad- 



vected with the flow of the explosion which record trajectories of temperature, 
energy, and density. Both this tracer method and the expanding grid add La- 
grangian components to our Eulerian formulation of the problem. 



4. Results 



Based on the numerical techniques described above, thermonuclear supernova 
explosion simulations in multiple spatial dimensions have been performed. While 
pure turbulent deflagrations are able to explode the WD, the resu lting event is 



predic ted to be on the faint end of the observed normal SNe la (jRopke et al 
20023) ■ 



The delayed detonation scenario, in contrast, covers the range of normal to 

brig ht SNe la when the ignition of the deflagration is treated as a s tochastic pro- 

cess ( Ropke fc Niemeverl 2007 [ Mazzali et al.ll2007l ) . In cases of few ( Garcfa-Senz fc Bravo 



20051-lL ivne et al.ll2005l: IRopke et al.l|2006bl) or asymmetrically distributed ign i 



.iVL 

tion sparks (e.g. ICalder et al.ll2004l : IRopke et al.1 l20"07bl : iTownslev et~aD 120071 ) . 



the deflagration is weak. It burns comparatively little material and hence the 
production of 56 Ni and the energy release are low. Therefore, when the deto- 
nation triggers, it finds abundant unburnt material at high densities which it 
converts to additional 56 Ni. The resulting event is bright and energetic and 
the ejecta structure is set mainly by the detonation stage. Somewhat counter- 
intuitively, a vigorous ignition in many sparks distributed around the WD's 
center and a subsequent strong deflagration leads to an overall faint explosion. 
Here, the detonation finds an almost completely burnt core and outer layers 
that have been diluted by expansion. It therefore burns them predominantly 
to IMEs and does not significantly contribute to the 56 Ni production. Hence, 
the structure of the iron group material at the center of the ejecta is dominated 
by the large-scale buoyancy instabilities from the deflagration. The outer IME 
layers, however, are produced by the detonation and therefore smooth. 
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A thorough exploration of the deflagration ignition configuration as a pa- 
rameter of the delayed detonation scenar io was performed in a set of two-di- 
mensional simulations ( Kasen et al. 20091 ). Here, spherical ignition kernels of 
radius 6 km were placed close to the center of the star - in some models isotrop- 
ically distributed around it and in others in a solid angle with an opening of 
less than 360°. In radius, a Gaussian distribution with a standard deviation of 
150 km or 75 km was chosen. The number of ig nition kernels rang ed from 15 to 
150 (see the Supplementary Online Material of lKasen et al.ll2009t ). In addition 



to the ignition spark distribution, the DDT criterion was varied. The detona- 
tion was triggered in the distributed burning regime if the ratio of the turbulent 
velocity fluctuation to the flame speed exceeded a certain threshold. Prom the 
results of these simulations, synthetic observables were computed by means of 
radiative transfer calculations. These found good agreement of the models with 
observations both in color light curves and in spectra and their evolution. Fur- 
thermore, from these results the peak brightness and the decline rate of the 
.B-band light curve could be determined for the individual mode ls. These were 
found to follow the observational relation of Phillips et al. ( 19991 ). 



5. Conclusions 

Simulating the turbulent combustion in thermonuclear supernovae is a numeri- 
cally challenging task. Specific techniques are required to correctly represent the 
thin flame and its interaction with turbulence on a wide range of spatial scales on 
a computational grid that comprises the full exploding WD. Here, a combination 
of a level-set based flame representation and subgrid-scale turbulence modeling 
were discussed as a solution to this problem. The implementation of this ap- 
proach allows to simulate thermonuclear supernova explosions. A pure turbulent 
deflagration is found to be inconsistent with the properties of "normal" SNe la. 
These objects require a detonation stage following an initial deflagration. The 
detonation wave is again represented numerically with the level-set technique. 
Simulations of the delayed detonation scenario lead to synthetic observables in 
good agreement with the observations. If the ignition of the deflagration is mod- 
eled as a stochastic process, the diversity of luminosities found in normal SNe la 
is reproduced. Moreover, the correlation between peak luminosity and decline 
rate of the 5-band light curve is reproduced. This is the first time that mul- 
tidimensional hydrodynamical explosion models predict this relation. Since it 
forms the basis of the calibration of SNe la as cosmological distance indicators, a 
theoretical understanding of the correlation is required to improve the precision 
of SN la cosmology. Our model provides the first step in this direction. 
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